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A coarse grained model for flexible polymers end-grafted to repulsive spherical 

nanoparticles is studied for various chain lengths and grafting densities under good 

solvent conditions, by Molecular Dynamics methods and density functional theory. 
O ' 

With increasing chain length the monomer density profile exhibits a crossover to the 

star polymer limit. The distribution of polymer ends and the linear dimensions of 

individual polymer chains are obtained, while the inhomogeneous stretching of the 

m ! 

CN) ■ chains is characterized by the local persistence lengths. The results on the structure 



factor of both single chain and full spherical brush as well as the range of applicability 
of the different theoretical tools are presented. Eventually an outlook on experiments 
is given. 
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I. INTRODUCTION 

Polymer chains densely grafted by a special endgroup to a substrate surface on which they 
otherwise do not adsorb, stick away from the substrate, forming a polymer brush.-"- Poly- 
mer brushes play a role of great scientific and technological relevance, even beyond colloidal 
stabilization.-^ The applications range from lubrication, tuning of adhesion and wetting 
properties, to biotechnology including protein purification, enzyme immobilization, virus 
capturing, improvement of biocompatibility of drugs, etc.— >ii - — The rich variety of practical 
uses in our daily life depends on the specific architecture and intramolecular interactions 
via several parameters, such as molecular weight of the chains, surface density, matching 
between the properties of the solvent and the monomeric units of the polymer. Due to the 
interplay between excluded volume interactions among the chains and the entropic repul- 
sion of the grafting substrate, polymer brushes acquire nontrivial structures where a polymer 
conformation is characterized by multiple length scales.— In this context the study of flat 
planar grafting surfaces has found longstanding interest, from the point of view of statistical 
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29l ). For many applica- 



mechanics of the macromolecular configurations (see e.g. Refs. 
tions (building blocks of nanocomposites^— , surface modification of biomolecules^^ etc.) 
the grafting surfaces are (approximately) spherical, with a radius comparable to the linear 
dimensions of the grafted polymers. However the effect of a spherical geometry of the sub- 
strate on the brush structure has been considered mostly for the extremely high curvature 
regime. It is indeed well known that when the radius of the central nanoparticle gets smaller 
than the typical size of the free macromolecule in solution, a crossover to the behavior of star 
polymers^— is expected. The intermediate curvature regime, where the size of the core is 
comparable to the full molecule extension, has received somewhat less attention.— >^r— An- 
other debated aspect is related to the description of micelles^ in terms of spherical polymer 
brushes. Small spherical nanoparticles with many flexible chains grafted to their surface are 
similar to spherical micelles, formed from asymmetric block copolymers AiBi_i with / <<1 
in a selective solvent (which is bad for the A-block). The A-monomers form a dense core, 



and the A-B interface can be idealized as a sphere, from which the B-blocks stick out^~— . 
In these models both the shape fluctuations of the A-rich core and the nonzero width of 
the A-B surface are ignored. Thus, an improved understanding of spherical brushes may 
be useful both for various applications of these nanoparticles, and to elucidate properties of 
related systems such as spherical micelles. 

In the present work, we study an idealized coarse-grained model for spherical polymer 
brushes both by Molecular Dynamics (MD) computer simulations and by calculations using 
density functional theory— We have / flexible polymer chains containing N effective beads, 
described by the standard Kremer-Grest^ 1 ^ model. These chains are grafted to sites regu- 
larly distributed on the surface of spheres such that the values of <r=0.068, 0.118 and 0.185 
of the grafting density are realized. The motivation for choosing these values originates from 
the fact that the same model and parameters have already been studied for planar brushes. 61 
The main purpose of the present work is to give a detailed quantitative description of the 
spherical brush behavior when the radius of gyration is comparable with the size of the core. 
At the same time, by changing the length of the polymer chains we explore the crossover 
between the flat brush regime (small curvature, short chains) and the star polymer regime 
(big curvature, long chains). More details on the model and simulation technique will be 
given in the following section ((IT])- In Sec JIHt we describe our numerical results, while Sec JIVI 
presents a discussion where comparison of MD to density functional calculations as well as 
simulation results for flat planar brushes are given. Eventually, referring to simulations and 
experimental literature we discuss the issue of which features of our results are universal or 
model-dependent. Sec JVl presents our conclusions. 

II. SOME COMMENT ON THE MODEL AND SIMULATION 
TECHNIQUES 

In order to realize the grafting densities <r=0.068, 0.118, 0.185 61 on spheres with a hard- 
core radius R c that is of comparable size to the gyration radius R g of the chain, we decided 



to choose respectively /=42 and R c =7, /=92 and R c =7.9, /=162 and _R C =8.35. With this 
choice of parameters R g can be conveniently calculated in our simulation. Here all the 
lengths are measured in units of the Lennard- Jones diameter o^j—l of the beads that form 
the chains. 

For planar brushes it is known^ that it makes little difference whether the grafting sites 
are regularly arranged at the grafting surface or at random, according to the chosen value of 
a. However to average out sample to sample fluctuations, the random distribution requires 
averaging over a large number of realization of grafting sites. To avoid unnecessary huge 
computation efforts, we decided to consider a regular distribution of the grafted monomers 
only. The method to find the appropriate positions of these grafting sites is a geodesic 
subdivision, i.e. a repeated subdivision of triangles. According to this approach one may 
begin with an icosahedron inscribed in a unit sphere, find the midpoint of each edge and 
renormalize its coordinates, pushing it out to obtain a new vertex lying on the unit sphere. 
This divides each original triangle into four new smaller triangles. It is possible to repeat 
this process and thus increase the number of vertices. One should note, however, that the 
lengths of the edges connecting the vertices are not all exactly the same, because of the 
projection to the unitary sphere. When the desired number of grafting sites is reached, the 
corresponding value of R c is readily calculated from a. 

The first monomer of each chains is rigidly fixed to the core. All the monomers interact 
via a truncated and shifted Lennard- Jones (LJ) potential: 
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for r < 2 1 / 6 a LJ 

(1) 
for r > 2 1 /Vj 



where r is the distance between the beads, while ctlj(=1) and £lj(=1) set the scales for 
length and energy respectively. The chains connectivity then is modelled by the finitely 
extendable nonlinear elastic (FENE) potential.— 41a2&i2Sl 

{ %ENE ( RsL) 2 l n \l_ (^ 2 

Vfene(^) 



for r < Rn 

(2) 
oo for r > R 



Here &fene = SO^lj and Ro = l.bau as usual.— di^M For this choice the total potential 
between bonded monomers (V=V^-|-Vfene) has a minimum at r~0.97o"Lj. EqfTJ mimics 
good solvent conditions (the solvent molecules are not taken explicitly into account). In our 
analysis we considered several chain lengths: iV=20,40,60 and 80 (the values do not take 
into account the grafted monomers). 

Standard MD methods using the velocity Verlet algorithm— were used, following previous 
works where multi-arm star polymers were simulated.— &L&L In order to stabilize the temper- 



ature of the system a Langevin thermostat as described in Refs. 



40 



.41 



.59 



.60 



is implemented 



and the chosen temperature is /c#T = 1.2. Assigning to the monomer mass m the value 



m = l, the characteristic time r = Jma^/e^ is then unity too. The integration time step 
was covering the range [10~ 3 r, 10 _2 r] with a total number of 1-2 xlO 6 timesteps used for 
equilibration and several independent runs of 3-5 xlO 7 timesteps each to gather statistics. 
FigJT] represents an illustrative example for the configurations thus generated. In the picture 
we shows a snapshot of a spherical polymer brush with / = 42, <r=0.068, and N=80. 

III. SIMULATION RESULTS 

We start our study analysing the conformational properties of the molecule. Figj2] shows 
the log-log plots for the radial monomer density profiles. Note that the grafted monomer 
at r=R c shows up as a delta-function spike. The pronounced structure close to the core is 
evidenced even further. The first mobile monomer (at about r=R c +ri lon d 1 r^nd ~0.97) is 
tightly bound to the grafting surface and hence shows up as a rather sharp peak. A further 
peak for the second monomer is also clearly seen, while thereafter (for N> 40) already 
the characteristic power law decay p(r) oc r" 4 / 3 expected for star polymers with very large 
arms 35 se £ s j n Ap ar t from a shift of the abscissa scale the data closely resemble simulation 
results for star polymers 4 ^ 1 ^. 

The pronounced peaks near r=R c + rb on d, are reminiscent of the layering of particles in 
fluids near hard walls, and this feature of our results is similar to the data for the density of 




FIG. 1. Snapshot of a spherical brush with /=42, <r=0.068 and iV=80. The additional grafted 
monomer is shown in light gray (yellow online). On the top (left side): the distribution of the 
grafted monomers on the sphere is shown. The free end monomers are highlighted as grey spheres 
(light blue online). 



polymer brushes grafted to flat planar walls.— i2&£I Of course in the region where the layering 
has died out, for flat polymer brushes one observes a decay of p(r) roughly compatible with 
a parabolic decay^i rather than the Daoud-Cotton power law observed here. For low enough 
N and/or big enough radius of the internal core, one should expect a crossover to the flat 
brush limit. As we will show in detail in Sec JIVI the results we obtained for iV<40 evidenced 
a parabolic profile. 

Finally, when the radial distance is r ~ R c + R gr (N) we observe a crossover of p(r) to an 
exponential decay (R gr (N) here is the gyration radius component of a chain in the radial 
direction.) 

Figj3] shows the corresponding density profile of the free chain ends. Unlike the Daoud- 
Cotton^ picture, chain ends do not only occur in the outer region of the spherical polymer 
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FIG. 2. Log-log plot for the monomer density p(r) at distance r from the center of the sphere of 
a spherical brush with /=42 (a), 92 (b) and 162 (c). The corresponding sphere radii and grafting 
densities are i? c =7.0, cr=0.068 (a), i? c =7.9, cr=0.118 (b) and i? c =8.35, cr=0.185 (c) respectively. 
Four choices of chain length, iV=20, 40, 60,80, are plotted in each figure, as indicated. The straight 
line illustrates the Daoud-Cotton 35 power law proposed for star polymers p(r) oc r~ 4 '' 3 '. 

brush, but rather throughout the whole brush. We point out here that for star polymers, 
scaling theories and renormalization group predictions 3 ^ 1 ^ imply 



p e (r) oc r &[}) ,r bond « r « R gr (N), 



(3) 



The exponent 9 can be expressed in terms of the exponents 7 and v of linear polymers 
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FIG. 3. Linear-linear plots of the free end monomer density profiles p e (j"), versus the radial distance 
r for the case /=42(a), 92(b), 162(c) and several chain lengths, as indicated. Insets (a), (b), (c) 
show the same data on log-log plots. 

and the partition function exponent 7(/) of / arm star polymers as 64 



0(/) = [7-7(/ + l) + 7(/)-l]/"- 



(4) 



Since for large / one expects 'j(f) oc — / 3//2 , Eq. 0] would imply a very large exponent, 
9(f) oc / 1//2 /V, which practically means that the chain ends are excluded from a substantial 
fraction of the internal corona of the star polymer. In contrast, for a brush on a planar 
grafting surface one finds from the self-consistent field theory in the strong stretching limit 21 



p e (z) oc z,r bond « z « Rgz(N), 



(5) 



where z is the coordinate perpendicular to the grafting surface. Unfortunately, R gr (N) for 
the available chain lengths is clearly not large enough to provide a significant regime of r 
where Eqj3] could be tested. The pronounced curvature of p e (r) on the log-log plots (inset 
in figure [3]) prevents us from making any strong statement on the exponent 9(f). In fact, 
on general grounds, one would expect a crossover from a linear behavior, when / is very 
large and N relatively short (such that R gr (N) << R c ) and the system locally resembles a 
polymer brush on a flat surface, to a very steep power law (Eqj3]) in the inverse limit, where 
R gr (N) » R c and the system resembles a star polymer. Our data, of course, can probe 
only a small intermediate part of this very extended crossover regime. Qualitatively, our 
results for p e (r) are similar to the results by Toral and Chakrabarti 46 for a different model 
(the so-called pearl- necklace model). We recall that for star polymers in d = 2 dimensions 
where 9(f) is known exactly, since 39 7(/) = 17/16+(9/32)[/ — /(/ — l)/2] is known from 
conformal invariance methods, Monte Carlo results for star polymers even for small / did 
not confirm Eqs. [3|4l 64 Clearly, it would be desirable to vary both / and N over a much 
wider range, but this is still impossible with the present method. 

Another issue of interest is the local stretching of the chains. Since the end-to-end vector 
R e of the chains, by the symmetry of the problem, is oriented on average in the perpendicular 
direction, we monitored the local persistence length l p (k) of the chains, defined as^ 

l P (k)/r bond = (at ■ R e /\dt\ 2 ), at = r£ - rv^, k=l,...,N. (6) 

Here rt denotes the position of the fe'th monomeric unit along a chain, starting with k=0 
at the grafting site. While for free chains (under good solvent conditions) l p (k) is known to 
exhibit a broad maximum near k=N/2, which scales as^ 1 ^ l P ,max(k) oc iV 2l/_1 , we find here 
(FigHJ) that l Pj max(k) occurs near the grafting site, for k=l, and decreases monotonically 
with k, getting of order unity at the free chain end. We recall here that we work with a fully 
flexible model, lacking any intrinsic stiffness of "chemical origin". Moreover under good 
solvent conditions we can not have quantitatively precise information on the persistence 



length as calculated from Eq. [6] (see Ref. 



68l ). However, Eq. [6] still gives a meaningful 



14 
12 
10 



f=42 0=0.068 



8'^r 



■ ■ N=20 
aan=40 
xx N=60 
• •N=80 



xxvvv ■ '"WMH ... 

A A A A A A A A A A A A A A A A A A A A A A A*!** 1**** 

■ ■ ■•■■■■ ■ ■■■■■■ ■ ■ ■••■■■ ■ ■ *:. m:~:^ 



(a) 



0.2 



0.4 0.6 

k/N 




15 



<H 10' 



f=92a=0.118 




^P^sxsrsswxawxK 



>***, 



,'iC.w 

* "*-» 



(b) 



0.2 



0.4 0.6 

k/N 



0.8 



N=60 0=0.118 

*■■■* N=60a=0.15 

N=20a=0.118 

*■•■-* N=20 a=0.15 




(c) 



0.4 0.6 

k/N 



FIG. 4. Local persistence length (as a measure of intrinsic chain orientation along the end-to-end 
vector) l p (k), Eq.[(6)], plotted vs. the normalized monomer index k/N for the case /=42(a) and 
92(b). Case (c) compares for /=92 two grafting densities (<r=0.118 for R c =7.9 and <r=0.15 for 
R c =7, respectively) to demonstrate that the dependence of l p (k) on grafting density is very weak, 
at least in the range of parameters investigated in this paper. 

information on the relative orientation of bond vectors, correlated with the end-to-end vector 
orientation. Nevertheless Fig. H] demonstrates that the behavior is mostly controlled by the 
excluded volume interactions between the monomers of the single chain; the repulsion due 
to the interchain interactions does not yet lead to a significant increase in orientation of the 
chains. On the other hand the chain linear size grows moderately in the radial directions 
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FIG. 5. Log-log plot of mean square radius and end-to-end distance components in radial and 
tangential directions versus chain length. Here R gr =J(R'i z ), R g t=J (R"i x + Ri y )/2, R er =J (B% z ), 



R e t=J (Rex + -^ej/)/2 orienting the z axis from the center of mass of a chain perpendicular to the 
grafting surface of the sphere, and the x, y axes are oriented perpendicular to this z-axis in each 
chain configuration that is analyzed. Two choices of / are shown, /=42 and 92. 



as the grafting density is increased. In Fig. [5] we show as an example the log-log plot of 
mean square radius and end-to-end distance components in radial and tangential directions 
versus N, for cr=0.068 and o = 0.118. The tangential components decrease with increasing 
/, in agreement with the increase of the stretching of the chains with the grafting density. 
Of course the differences reduce as the number of monomers grows. Note, however, that all 
the chain linear dimensions still scale with chain length according to the standard excluded 
volume power law, R(xN v . Neither the scaling relation for star polymers 35 



(EEL) oc (Rl) ex JV 5 */ 1 -" 



(7) 
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nor for chains in planar brushes in the strong stretching limit 21 

{R%) « (Rl) °c iV 2 V/ 3 (8) 

describe our numerical results perfectly. Thus, the conclusion again is that there exists a 
very extended regime of gradual crossover from the planar brushes to the star polymer limit. 
The results we obtained by MD simulations probe only a small section of this regime, as 
noted above. 

It is worth pointing out that R gr and R g t are of the same order of magnitude; if one 
would take blob pictures like that of Daoud and Cotton^ literally, one would conclude each 
chain "lives" in an angular sector of extent 47r//, and hence for large / one would have 
Rgt « Rgr- Obviously, for entropic reasons it is more favorable that the different chains 
share these angular sectors, and thus the chains are also less stretched in radial directions 
than expected from the Daoud Cotton picture. 

As a last point of this section, we consider the structure factor describing the total 
scattering from the spherical brush, as well as the structure factor describing the spheri- 
cally averaged scattering intensity from single arms (Figj6]). While the latter resembles the 
structure factor of an isolated polymer under good solvent conditions crossing over from 
S(q) = N(l — q 2 < R 2 g > /3) at small q to S(q) oc q^ 1 ^ at larger g 65 , the structure fac- 
tor for the full spherical brush clearly exhibits a more interesting behavior. It resembles 
the scattering from a compact sphere at small q and exhibits a minimum, the location of 
which is given roughly by 1/R C ] the maximum at larger q is more and more washed out as 
N increases. These data clearly resemble quantitatively experimental scattering data from 
spherical micelles.— 1 ^ 

IV. COMPARISON TO RESULTS FROM DENSITY FUNCTIONAL 
CALCULATIONS AND TO OTHER RELATED WORKS 

In this section we mainly focus on three different aspects: the limit of applicability of the 
MD method, the density functional theory as a possible approach to extend the accessible 
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FIG. 6. Log-log plot of the structure factor of a spherical brush (empty circles and full circles) 
and of a single arms (full curves without symbols) plotted versus wavenumber q, in the case R c =7, 
<j=0.068, /=42, and R c =7.9, <r=0.118, /=92, for two chain lengths: iV=60(a) and iV=80(b). 
Broken straight lines indicate the power law S(q) oc q^ 1 or S(q) oc q^ 1 ' 1 ' = g- 1 / - 588 ^ respectively. 

range of parameters, the peculiarities of spherical brushes (crossover regime) with respect 
to star polymers and planar brushes. 

In view of the rather large computational efforts needed in order to generate well equi- 
librated MD data on spherical polymer brushes, with very good statistical accuracy and 
varying a number of parameters, it is worthwhile to ask whether or not one can obtain 
results of comparable quality using a "cheaper" semi-analytical method, such as density 
functional theory (DFT).-3S>S& 

Since the details of these calculations have been given in a different context elsewhere^ 



13 



i 1 — i — i — i — i 1 — r 



N=80 MD 
N=60 MD 
N=40 MD 
N=20 MD 




FIG. 7. Log-log plot of the monomer density p(r) versus r for /=42(a) 92(b), 162(c) comparing 
MD results (broken lines) with DFT results (full lines), for iV=20, 40, 60 and 80 (from left to 
right). In the picture the results for 7V=250, 500, 750 obtained via DFT are also shown (the three 
right most curves). 

we here describe the results only, but we emphasize that a perfect agreement can not be 
expected, since the DFT uses a slightly different model for the chain molecules, employing 
tangent hard spheres rather than FENE+LJ potential (Eqs. [TfiZj) to create chain connectivity. 
As a result, the chains molecules modelled by DFT are slightly different that those studied 
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FIG. 8. Log-log plot of the free end monomer density p e {r) versus r for /=42(a) 92(b), 162(c) 
comparing MD results (broken lines) with DFT results (full lines), for iV=20, 40, 60 and 80 (from 
left to right). 

by MD, leading to a slight mismatch between the results of both methods. Nevertheless one 
recognizes that the DFT reproduces the general trend of the MD results rather faithfully 
(Fig. IS]). Also the layering of p(r) close to the grafting surface is accurately reproduced 

(Fig. ED- 

In order to study the crossover of the average brush height H av from the flat brush limit 
to the star polymer limit in detail, we follow Wijmans and Zhulina 45 by defining H av from 
the second moment of the monomer density profile as follows 



H: 



. oo 
2 / J2 



2 / / „2 



H c 



r drp(r)(r — R c ) / r drp(r) 



R c 



(9) 
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FIG. 9. Log- log plot of the normalized brush height H av /Ho versus the normalized inverse radius 
Hq/R c , comparing MD results (open symbols) to DFT results (filled triangles, a = 0.185; crosses, 
a = 0.118; asterisks, a = 0.068). The theoretical asymptotic slope predicted from SCFT (-2/5) is 
shown as a broken straight line. The DFT results converge to a power law with slope of minus one 
third (or -2/6), (full straight line). Note that H av was computed from Eq. [9] while the corresponding 
height Hq of a planar brush, using the same values of chain length iV and grafting density a, we 
calculated from Eq. [TUJ 



The limiting height for a planar brush, obtained from Eq. [9] as R c — > oo, is denoted as 



Ho = / p[z)z 2 dzj \ dz p(z) 
Jo Jo 



(10) 



From the SCFT treatment of Wijmans and Zhulina^ we conclude that in the limit R c — > oo, 
iV — > oo (these limits need to be taken together such that < H Q /R C < oo) a scaling property 
holds, where /(^ a ) is some scaling function, 



H av / Hq = f(H /R c 



(11) 



irrespective of R c and N. The scaling is also irrespective of the grafting density a, provided 
it is chosen such that the monomer density in the brush stays in the semidilute regime: the 
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crossovers to either the mushroom regime or to dense brushes where p(r) attains melt-like 
densities near the grafting surface are not covered. 

As it is well-known, SCFT in the flat brush regime predicts that£~— 

H = a l/3 N (12) 

On the other hand, the star polymer scaling limit 35 implies that H av oc A 3 / 5 . As a conse- 
quence, the scaling function f(X) for large X must behave as /(A) oc A~ 2 / 5 , and hence one 
predicts 3 -^^ 5 -^ 

H av oc H 3/5 R 2 J 5 oc a^N^Rl'K (13) 

FiglH] presents a plot of our MD data, using such a scaling representation, and compares 
them to the corresponding DFT results. The latter method has the distinctive advantage 
that it can be applied to much longer chains. While (for our choices of R c ) the MD data do 
not extend beyond H /R c « 2.5, the DFT data can easily be taken to H /R c w 20, using 
chain lengths up to A=750, which would be inaccessible for our MD study. As expected 
from the comparison of the profiles, DFT and MD results are close to each other, and exhibit 
a similar trend. However, it is also clear that our data do not fall in a regime where the 
scaling described by Eq. [TT] is valid: the results for the three grafting densities used in our 
MD work do not superimpose on a single "master curve", which then would yield an estimate 
of a piece of the scaling function /(A). It is also clear that the MD data did not reach the 
regime where Eq. [13] is valid. The DFT results, on the other hand, do reach an asymptotic 
power law but it is inconsistent with Eq. [131 rather implying that f(X) oc A~ 1//3 . We lack 
an explanation as to why DFT predicts a star polymer scaling R s t ar oc JV 2//3 (following from 
f(X) oc A^ 1 / 3 ) rather than R star oc A 3 / 5 . 

These discrepancies are also apparent when we study the end-to-end distance of the 
chains, predicted by DFT for both flat and spherical brushes (Fig. fTOj) : while for flat brushes 
DFT predicts correctly R e oc A, for spherical brushes one finds R e oc A 2 / 3 rather than the 
scaling result, R e oc A 3 / 5 . Thus DFT is less useful for spherical brushes than for flat brushes, 
overestimating the chain stretching in spherical brushes for large A systematically. 
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FIG. 10. Log- log plot of the end-to-end distance of polymers grafted in flat and spherical brushes 
as obtained from MD (full symbols) and DFT (open symbols); straight lines indicate exponents 1 
(broken straight line), 2/3 (dash-dotted straight line) and 3/5 (full line), respectively. Inset: Log- 
log plot of the end-to-end rescaled distance Re/a 1 ' 5 versus the number of monomers per chain, for 
iV<80. 

It is also of interest to directly compare the density profiles as obtained from MD for 
flat brushes and for spherical brushes (Fig. [TT]) . As expected, for N = 20 the differences 
between the profiles for flat and spherical brushes are rather minor, since for the chain values 
of R c then H av << R c for all these grafting densities. On the other hand, for N = 60 the 
curvature of the grafting surface has already a clearly recognizable effect on the profile. 

For completeness we plot the radial and tangential components of the radius of gyration 
and end-to-end distance of the single polymer grafted to the flat brush. The figure shows the 
different scaling trends with respect to the spherical brush cases (see Fig. [12]for comparison). 

Finally we mention that in the recent experiment of Dukes et al. 32 a scaling H av ex 
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FIG. 11. Monomer density profiles p(r) for iV=20 and iV=60 plotted versus distance r from the 
grafting surface for both spherical brushes and flat brushes, as indicated, including three grafting 
densities: a = 0.068 (upper panel), a = 0.118 (middle panel), a = 0.185 (lower panel). All 
the profiles are normalized to unit area /q°° dz p(z)=l (flat brushes) or JS° dr r 2 p(r)=l (spherical 
brushes), respectively. 

TV 4 / 5 was observed, eventually crossing over to H av oc N 3 ^ 5 at larger N. The authors 
emphasize, however, that this behavior should not be explained in terms of Eq. [11] since their 
brushes are in a regime where near the grafting surface the semidilute concentration regime 
is distinctly exceeded, and one must build the description on a theory for concentrated rather 
than semidilute polymer brushes. Nevertheless they find that the approximate prediction 
of Wijmans and Zhulina^ for the crossover grafting function f(X) roughly accounts for 
their data, as well as for previous data of other authors^^i. However, it is clear that a 
fully satisfactory theoretical description of these data does not yet exist. Nevertheless it is 
encouraging that experiments have become feasible for conditions that are not so different 
from the conditions studied in the present work. 
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FIG. 12. Log- log plot of mean square radius and end-to-end distance components in radial and 
tangential directions versus chain length. Here R gr =J(R g ) z ), R gt =, / (Rg X + R% y )/2, R er =J (R% z ), 



R e t=J(Rg X + Rl y )/2 orienting the z axis perpendicular to the flat grafting surface, and the x,y 
axes lie on the plan parallel to the surface in each chain configuration that is analyzed. Two choices 
of / are shown, / =42 and 92. 

V. CONCLUSIONS 

In this work we have presented a detailed and quantitative analysis of conformation and 
structural properties of spherical brushes. 

In the first instance we focused on molecules where the core has a radius comparable to 
the polymer linear dimension. In particular we have carried out monomer resolved Molecular 
Dynamic simulations in order to study the monomer density profile and the distribution of 
the free end monomers for three values of the grafting density (namely cr=0. 068, 0.118, 0.185). 
In order to clarify the origin of the inhomogeneous stretching of the chains, and the devel- 
opment of multiple relevant length scales we calculated the local persistence length, the 
structure factor of single chain and full molecule as well as the radial and tangential compo- 
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nents of the radius of gyration and end-to-end distance of the single polymer. In the regime 
of parameters investigated the structure of the brush is mostly driven by the intrachain ex- 
cluded volume interactions. The chains are less stretched in radial directions than expected 
from the Daoud Cotton picture. The main conclusion here is that there exists an extended 
regime of gradual crossover between star polymer and flat brush behavior. Using the ratio 
of the brush height at a flat surface to the core radius as a variable for this crossover, we 
can state that this crossover extends over several decades in this variable, in order to cover 
both star polymer scaling and flat brush scaling relations. In addition, one has to reach 
the strong stretching limit of semidilute polymer brushes (avoiding both mushroom-like be- 
havior and too dense melt-like brushes) to observe this simple universal crossover scaling 
behavior. Both by simulations and in typical experimental systems it is not easy to reach 
this asymptotic scaling regime, and hence the crossover scaling considerations are only a 
rough guidance to interpret the expected behavior. 

We also comment on the use of the present model to interpret scattering from spherical 
micelles. Indeed the chain lengths and total number of chains grafted to the spherical 
core particle are in a similar range as for typical micelles in experimental studies. While 
the scattering curves found here (jSJ) are quantitatively similar to experimental data, we do 
expect some differences due to the fact that in micelles there is a smooth interfacial profile 
between the monomers of type A in the collapsed core and the monomers of type B in the 
(swollen) corona (referring here to block copolymer micelles of composition AiB\_i with 
/ <<1 in a selective solvent, poor for A but good for B, as used in typical experiments). 
Therefore one cannot expect the pronounced layering in the density profile of the corona 
monomers, that we see here for radii r, near the radius R c of the hard particle to which our 
chains are grafted (JTJ). 

In the second instance we considered the density functional approach as a possible al- 
ternative method. The results obtained via DFT are in a very good agreement with our 
simulation results and open up the possibility to extend the range of parameters investigated. 
Nevertheless the theory does not reproduce exactly the star polymer scaling. Moreover com- 
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paring the results obtained on flat brushes via MD simulations and DFT we observed that 
DFT is more precise in describing the scaling of flat than spherical brushes, overestimating 
the chain stretching as the radius of curvature of the surface increases. 

We finally compare the MD data for flat and spherical brushes in order to emphasize the 
peculiarity of the regime characterized by a grafting surface with a radius comparable to the 
polymer linear dimension. 

This crossover regime of the surface curvature radius (between flat interfaces or star-like 
polymers) is of great interest for many practical purpose closely related to our daily life. 
One aspect we are interested in is related to the investigation of the adsorption of proteins 
to polymeric spherical brushes, e.g. how adsorbed or chemically-attached polymers affect 
the interaction between nanoparticles and globular proteins. A full understanding of the 
structural properties of the polymeric unit is a first necessary steps in a such direction. 
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